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Abstract 

Appfications to holographic theories have led to some recent interest in magnetic 
monopoles in four-dimensional Anti-de Sitter spacetime. This paper is concerned with 
a study of these monopoles, using both analytic and numerical methods. An approx- 
imation is introduced in which the fields of a charge monopole are explicitly given 
in terms of a degree N rational map. Within this approximation, it is shown that the 
minimal energy monopole of charge has the same symmetry as the minimal energy 
Skyrmion with baryon number in Minkowski spacetime. Beyond charge two the 
minimal energy monopole has only a discrete symmetry, which is often Platonic. The 
rational map approximation provides an upper bound on the monopole energy and 
may be viewed as a smooth non-abelian refinement of the magnetic bag approxima- 
tion, to which it reverts under some additional approximations. The analytic results 
are supported by numerical solutions obtained from simulations of the non-abelian field 
theory. A similar analysis is performed on the monopole wall that emerges in the large 
A'^ limit, to reveal a hexagonal lattice as the minimal energy architecture. 
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1 Introduction 



The AdS/CFT correspondence allows the investigation of strongly coupled theories by 
studying classical solutions in the bulk. Recently, it has been argued that interesting phe- 
nomena may result, including spontaneous breaking of translational symmetry, if the bulk 
Anti-de Sitter (AdS) spacetime contains non-abelian magnetic monopoles [1]. This pro- 
vides motivation for a detailed study of SU{2) magnetic monopoles in four- dimensional AdS 
spacetime. 

In addition to holographic applications, there are other reasons to consider monopoles in 
AdS. In Minkowski spacetime, when the Higgs field is massless, there is a 4A^-dimensional 
moduli space A^at of static charge BPS monopoles. In contrast, the attraction in AdS 
spacetime should produce a unique minimal energy monopole for each charge N (up to the 
obvious action of spatial rotations). It is of interest to determine the structure and symmetry 
of the A^-monopole, which will map to a unique point in the moduli space A^at as the AdS 
curvature tends to zero. 

It has been known for some time that there are many similarities between BPS mag- 
netic monopoles and Skyrmions (for a review see [2]). However, one important difference 
is the attractive force between Skyrmions, producing bound states, in contrast to the BPS 
monopole moduli space Mn, resulting from the absence of static forces between monopoles. 
It is therefore expected that studying monopoles in AdS will enhance the similarities with 
Skyrmions, since both will share the features of attractive forces and bound states. 

Motivated by previous work on Skyrmions [3], an approximation is introduced in which 
the non-abelian fields of an A^-monopole in AdS are written in terms of a degree N rational 
map between Riemann spheres. The Yang-Mills-Higgs energy functional leads to an energy 
functional on the space of rational maps, which is precisely the one found in the similar 
approach to Skyrmions. This implies that, within this approximation, the minimal energy 
monopole of charge N in AdS has the same rotational symmetry group as the minimal energy 
Skyrmion with baryon number and massless pions in Minkowski spacetime. 

In particular, this approach predicts that the 1-monopole is spherically symmetric and 
the 2-monopole is axially symmetric, but for charges greater than two the minimal energy 
monopole has only a discrete symmetry group, which is often Platonic. 

In Minkowski spacetime the existence of particular BPS monopole solutions with Platonic 
symmetries has been proved [H [S] using the Nahm transform [6J, but in the BPS case there 
is no energetic preference for these solutions, because of the democracy of the moduli space 
A^AT. Although the rational map approximation does not produce any exact solutions, 
beyond charge one, it does provide an upper bound on the A^-monopole energy and suggests 
that in moving from Minkowski spacetime to AdS, the breaking of the energy degeneracy 
of the moduli space Ai^ leaves particular symmetric monopoles as the minimal energy 
solutions. 

Monopoles with large charge have been considered previously using the magnetic bag 
approximation, both in Minkowski spacetime [7] and in AdS |T]. The magnetic bag approxi- 
mation assumes that inside a spherical bag both the Higgs field and the magnetic field vanish 
and outside the bag there is an abelian magnetic field. The rational map approximation may 
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be viewed as a smooth non-abelian refinement of tlie magnetic bag approximation, to wliicli 
it reverts under some additional approximations. 

Simulations of the full non-abelian field theory are performed to provide numerical so- 
lutions for monopoles in AdS, with charges from 1 to 17. These numerical results provide 
support for both the approximate analytic treatment using rational maps, and the magnetic 
bag approximation. The numerical results also reveal interesting information about the zeros 
of the Higgs field. For most of the solutions, there is a zero of the Higgs field with multi- 
plicity N located at the origin. However, for N = 3 and some further non-minimal solutions 
with higher charge, there are + 2 zeros of the Higgs field, with one of the zeros having a 
negative multiplicity. 

By zooming to the Poincare patch in the large N limit, the study of magnetic monopoles 
translates into the investigation of monopole walls in AdS. A monopole wall [8] is a novel 
domain wall in which the magnetic field along a line perpendicular to the wall tends to zero 
on one side of the wall and to a non-zero constant on the other side of the wall. The fields 
are periodic in the directions parallel to the wall, with a non-trivial spatial variation of the 
energy density and magnetic field. Monopole walls have infinite energy per unit area in 
Minkowski spacetime, but in AdS this is finite and they can be studied using a magnetic 
bag style approximation [IJ. 

The magnetic bag approximation is too crude to reveal any information concerning the 
spatial distribution of the fields parallel to the wall, but a variant of the rational map 
approximation, involving elliptic functions, is refined enough for this purpose. It suggests 
that the minimal energy per unit area is obtained from a monopole wall with a hexagonal 
lattice. Numerical simulations are performed that support this conclusion, which is again in 
agreement with the Skyrme model, where a similar hexagonal lattice exists [5]. 

Previous studies of monopoles in AdS have been restricted to the spherically symmetric 
single monopole [lOlITT], and axially symmetric monopoles with charges two and three |12] . 
The computations in p.2j reveal that the energy of the charge two monopole is greater 
than twice the energy of the single monopole, but this result was incorrectly interpreted 
to conclude that monopoles repel and that the axially symmetric 2-monopole is unstable. 
The source of the confusion is that in Minkowski spacetime the following argument can be 
applied: if the energy of an A^-soliton is greater than N times the energy of a single soliton, 
then the energy of the A^-soliton can be reduced by infinitely separating the constituents. 
However, in AdS this argument fails because the energy of a single soliton increases as it 
approaches the boundary of AdS. The energy of the single soliton, which is used in the 
false comparison, is the single soliton energy only when the soliton is located at the origin 
of AdS. At first glance this may appear confusing, given the large isometry group, 0(3,2) 
of AdS, but there is an important subtlety at work here. To consider static solutions, and 
their associated energies, requires the selection of a time-like Killing vector. This breaks the 
isometrics of AdS and introduces a preferred point in space, which is denoted the origin. 
Solitons are attracted towards the origin, which is consistent with the interpretation of AdS 
as gravitational attraction. 
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2 Monopoles, rational maps and magnetic bags 



In terms of sausage coordinates, the metric of four-dimensional AdS spacetime may be 
written as 

ds' = - ([^) ' dt' + ^^i^ l^dp' + p\de' + sin^ e V)) , (2.1) 

where the range of the radial coordinate is < p < 1 and L is the AdS radius, related 
to the cosmological constant via A = — S/L^. The above choice of coordinates for AdS is 
motivated by later numerical investigations, as it is more efficient to have a finite range for 
the coordinates. 

The Yang-Mills-Higgs action is 

S = JItt (^F^^F^^ + Id,<!>D^<^^ d'x, (2.2) 

where the gauge potential and Higgs field $ are both sM(2)-valued. The massless Higgs 
field transforms in the adjoint representation of SU (2) and is subject to the boundary con- 
dition |$p = -iTr$2 -)■ 1 as p-)- 1. 

At a fixed time, the Higgs field on the spatial boundary is a map between two spheres, 
$|p=i : "S*^ S"^, with winding number G 'n'2{S'^) = Z. This winding number is equal to 
the magnetic charge of the monopole, in units of 27r, and may also be written as the integral 

^ = ^ / £^Jk^TiFJkDi^) d^X. (2.3) 



Static monopoles are critical points of the static energy associated with the action (12. 2p . 
namely 

E = -j^Ti (^\f,,F'^ + ^DiW<!>^ ^ d'x. (2.4) 

The aim is to find the global minimum of this energy, within each topological sector given 
by the positive integer A^. Note that the above normalizations are such that in Minkowski 
spacetime each point of the BPS moduli space A^at is associated with a monopole solution 
with energy 2'7rA^. 

It is perhaps worth pointing out that although constant time slices of the metric ( 12. ip 
give the Poincare ball model of hyperbolic space, the static energy (12. 4p . and hence the static 
field equation, is not the much studied system for monopoles in hyperbolic space [13]. This 
is because of the extra factor y/—gtt = (1 + P^)/(l ~ P^) hidden in the term y/—g from the 
warp factor in the metric (12. ip . 

In Minkowski spacetime there is a diffeomorphism, preserving the action of spatial ro- 
tations, between the BPS moduli space Ai^ and a certain equivalence class of degree A^ 
rational maps between Riemann spheres [T4]. The rational map appears as scattering data 
along a radial half-line of an operator constructed from the monopole fields. As the inverse 
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scattering problem is not tractable for > 1, there is no explicit formula for the monopole 
fields in terms of the rational map. 

Motivated by similarities between monopoles and Skyrmions, an explicit but approximate 
description of Skyrmions was introduced in terms of rational maps |3]. Continuing the 
theme of an interplay between monopoles and Skyrmions, the rational map approximation 
for Skyrmions will now be adapted back to the monopole context, to provide approximate 
fields for monopoles in AdS. 

A rational map between Riemann spheres, R{z), is simply a ratio of two polynomials in 
a complex variable z. The two polynomials are required to have no common roots, and the 
degree of the rational map is the largest of the degrees of the two polynomials. To make 
the connection to monopoles (or Skyrmions) the Riemann sphere coordinate z is related 
to the angular space coordinates 9 and ip via standard stereographic projection, that is 
z = e'^tan{9/2). 

The Riemann sphere coordinate, R, on the target space of the rational map, is related 
to a three-component unit vector n by inverse stereographic projection, 

^ = TTT^(^ + R,t{R-R),l- (2.5) 
1 + \R\ 

The approximate monopole fields in AdS are taken to be 

^ = iHn-T, Aj = ^{l-K){nxdjn)-T, (2.6) 

where H{p) and K{p) are real radial profile functions. Regularity at the origin imposes the 
boundary conditions H{0) = and K{0) = 1. The conditions at the spatial boundary of 
AdS follow from the requirement that as p — )■ 1 then |$| — )■ 1 and Dj^ — )■ 0. This provides 
the boundary conditions H{1) = 1 and K{1) = 0. 

A more general form for the gauge potential, using all the symmetric tensors that can be 
computed from n, appears to be 

Aj = '-{1 - K){n X djiY) ■T+'-Pdjn-T+ '-Qx^n ■ r, (2.7) 

with additional profile functions -P(p) and Q{p). However, Q corresponds to an abelian gauge 
potential in the reduced one- dimensional radial theory, and therefore may be set to zero by a 
gauge transformation. Furthermore, there is a global U{1) symmetry that rotates the fields 
K and P, and this can be used to set P to zero. 

As $|p=i = m-T, it is obvious that the monopole charge N, which is the winding number 
of the Higgs field on the boundary two-sphere, is equal to the degree of the rational map 
R{z) that determines n via (12. 5p . It will be useful later to note the integral expression for 
the degree 
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where the final factor in (12.81) is simply the standard area element on the two-sphere. 



5 



Substituting the approximate fields (12. 6p into the energy (12.41) yields the expression 



i?™* = 4vrL {-(^^^2ip'H" + 2NH'K') + ^ (^2NK" + y,{K' - 1)^) } (1 + p') dp, 

^ ^ (2.9) 

where X is defined as the following energy functional on the space of rational maps 
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The functional (I2.10p is precisely the same quantity that appears in the rational map 
approximation of Skyrmions [3], and in both situations this is the only contribution to the 
energy that distinguishes between rational maps of the same degree A^. This proves that, 
within the rational map approximation, the minimal energy monopole of charge N in AdS 
has the same rotational symmetries as the minimal energy Skyrmion with baryon number 

and massless pions in Minkowski spacetime. 

Previous results on rational maps, obtained in the context of Skyrmions, will therefore 
be of use later in this paper. For example, an inequality that follows simply from (12.81) 
and (I2.10p is X > A^^ [3j. Furthermore, for a large range of A^, including all A^ < 40, 
the rational maps that minimize X have been computed numerically [3l [151 [IS] ) cind their 
symmetries identified. Recall that a rational map is symmetric under a group G G 5*0(3) 
if a spatial rotation g & G, which acts on the Riemann sphere coordinate z as an SU{2) 
Mobius transformation, can be compensated by an SU (2) Mobius transformation acting on 
the target Riemann sphere coordinate R. A rotation on the target sphere corresponds to a 
gauge transformation, and hence a symmetry, of the monopole. 

The monopole fields (12. 6p are consistent with the static field equations only if the rational 
map is spherically symmetric. The only spherically symmetric rational map is the (unique 
up to the action of SU{2) Mobius transformations) degree one map R = z. In this case 
X = A^ = 1 and the rational map approximation is exact, with the expression (12.90 for E^°'^ 
reproducing the true monopole energy. In this case the ordinary differential equations for the 
profile functions H and K, that follow from the variation of (12. 9p . agree with those appearing 
in the previous investigations [HI, llj of the single monopole in AdS. 

A numerical construction of the profile functions allows a computation of the 1-monopole 
energy as a function of the AdS radius L. The results reveal that for L > ^ the energy is well- 
approximated by the formula E ^ 27r(l + ^L~^). In particular, the energy is = 27r x 1.396 
for L = 1. In all the numerical computations presented in this paper the value L = 1 will 
be chosen as a generic radius. The above result shows that for this radius the energy of 
the single monopole is increased by around 40% from the fiat space limit, which should be 
sufficient to observe the phenomena that arise due to the curvature of AdS. In the analytic 
approximations the dependence on L will be retained, allowing the qualitative behaviour 
with L to be deduced. 

For each charge A^, using the minimal value of X in (12. 9 p and numerically computing the 
profile functions H and K, produces an energy E^"''^ which is an upper bound on the true 
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minimal energy of the A^-monopole. The energies E^"-^ (in units of 2tt) are plotted as the 
circles in Figure [T] for 1 < < 17. 




N 

Figure 1: The energy (in units of 2tt) for monopoles with charge N < 17. Circles denote 
the values obtained from the rational map approximation and squares are the results of 
field theory simulations. The solid curve is N^/"^ and the dashed curve is the magnetic bag 
approximation. 



The rational map approximation may be related to the magnetic bag approximation 
O [1] by making a series of further simplifications, as described below. The magnetic bag 
approximation assumes that the fields take different forms inside and outside a bag of radius 
p = a. Inside the bag, that is p e [0, a) the profile functions are taken to be if = and 
= 1, so that the energy density vanishes. For p > a, which is the outside of the bag, the 
magnetic field is taken to be abelian by setting K = 0. With these simplifications the energy 

Ei:l = 4nL [' { /^'l^ + -I^]{l + p^)dp, (2.11) 

where the contribution from the surface of the bag has been ignored, which means that E^^^ 
is no longer guaranteed to be an upper bound for the true monopole energy. The boundary 
condition on the surface of the bag is H{a) = 0, and the energy minimizing profile function 
is easily found to be 

H = ^^;^l^, (2.12) 
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where the function 

h{p) = - - p + AidjT^ p, (2.13) 
P 

has been introduced. The resuhing energy depends on the bag radius a via 

The magnetic bag approximation is expected to be accurate for large charge N ^ when 
the bag radius approaches the AdS boundary, that is, cr = 1 — e with < e <^ 1. To leading 
order in e 

^C = 2-(i + i). (2-15) 

and the energy is minimized when e = \/QL/X^/^, to give the value 




< = 2-^/^X3/^ -2vr^^ (2.16) 



where the final approximation uses the fact that, for a large range of A^, numerical results 
reveal [I6] that X 1.3A^^ ~ (|)^/^A^^. In fact, the analysis in section H] suggests the limiting 
behaviour X/iV^ — )■ 1.21 as — )■ oo, which is consistent with the numerical results for large 
N. 

The solid curve in Figured] is the function A^^/^, obtained from the approximation fl2.16p 
by setting L = 1. It can be seen that this approximation is a reasonable fit to the numerical 
data (circles) and confirms the superlinear growth of E^""^ with A^. As discussed in the 
introduction, in AdS a superlinear growth of the monopole energy with charge does not 
imply that the A^-monopole is unstable to fragmentation into individual monopoles. 

The original magnetic bag approximation |71 H] assumes that the angular distribution is 
spherical. In terms of the rational map approximation this corresponds to the assumption 
that the winding density 



1 + u 



1 + \R 



' (2.17) 



dR 



dz 



appearing in the expression (12. 8p for the charge, is constant. Of course, the only rational 
map for which this is true is the degree one map R = z, corresponding to the fact that there 
are no spherical monopoles with A^ > 1. However, ignoring this fact and assuming that the 
density f l2.17p is constant, then it must be equal to A^. In which case X = A^^, so the spherical 
assumption is equivalent to approximating X by its lower bound A^^. With this additional 
simplification the energy E"^^* in (I2.16P becomes the bag energy 



/2A^3 

^6a, = 27ry^ (2.18) 

first obtained by Bolognesi and Tong [1] . The dashed curve in Figure [1] is this magnetic bag 
energy. 
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The energy i?^'^* of the rational map approximation is an upper bound for the true 
monopole energy, and the bag energy Ehag is obtained by neglecting certain contributions to 
the energy and assuming an idealized spherical distribution. It might therefore be expected 
that the true monopole energy lies between these two approximations. In the following 
section numerical results will be presented that confirm this expectation, using simulations 
of the full non-abelian nonlinear field theory. 



3 Field theory simulations 

To perform numerical field theory simulations it is convenient to introduce Cartesian type 
coordinates x, defined inside the unit ball, so that the sausage metric becomes 

where = x • x < 1. 

The numerical monopole solutions are obtained by performing a simulated annealing 
energy minimization algorithm on the associated static energy 

The coordinates x are discretized on a regular lattice with lattice spacing dx = 0.02, and the 
simulation grid contains all the points of the lattice that satisfy x - x < 1. Spatial derivatives 
are approximated using a second order finite difference scheme and the energy is computed 
at points of the dual lattice. On the boundary of the grid the vacuum expectation value 
of the Higgs field is enforced |$| = 1, and $ is fixed at the sampled values of a continuum 
field with winding number N. Explicitly, on the boundary $ = in ■ r with n given in terms 
of a degree N rational map via (12.51) . In theory, all degree rational maps should provide 
equivalent boundary conditions, as the winding number of the Higgs field is the only gauge 
invariant quantity. However, on the lattice there is likely to be a small bias against rational 
maps with extreme angular derivatives. 

The final ingredient required for the numerical simulation is an initial condition. The 
only requirement on the initial fields is that the Higgs field on the boundary must take the 
form described above, namely $ = m • r, with n determined by a particular rational map. 
The simplest possibility is to take the initial fields to be 

$ = zpn ■ T, Aj = 0. (3.3) 

Note that in the continuum theory this initial field does not have finite energy, because Dj^ 
does not vanish on the boundary. However, this is not a problem on the lattice, and the 
gauge potential evolves during the simulation, and in particular at the boundary, to drive 
Dj^ towards zero at the boundary. A more sophisticated initial condition could be used. 
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Table 1: The symmetry group and energy (in units of 27r) for monopoles with charge < 17. 

for example by starting with the fields of the rational map approximation, but this is not 
necessary. 

An indication of the numerical accuracy of the algorithm can be obtained by computing 
the charge A^, of the final configuration, using the lattice version of (12. 3p . For the results 
presented in this section the numerical charge is integer- valued to an accuracy of within 0.1% 
for A^ < 4, with the error rising to 1% for A^ = 17, which is the largest value considered. 

The minimal energies obtained from field theory simulations are presented in Table [1] 
and plotted as the squares in Figure [H As expected, these energies lie between the values 
predicted by the rational map and magnetic bag approximations. The rational map approxi- 
mation is more accurate for small values of A^, but the magnetic bag approximation becomes 
increasingly accurate as A^ increases. 

To obtain these numerical results the rational map used in the initial (and boundary) 
condition is the X minimizing map for each value of A^. The symmetry of this map, which is 
also found to be the symmetry of the final numerical solution, is listed in Table [U Energy 
density isosurfaces are displayed in Figure [2] for a selection of monopoles with continuous or 
Platonic symmetries. For N > 2 the energy density is localized on the edges, and particularly 
the vertices, of a polyhedron. For the larger charges shown in Figure |2]the pattern emerges of 
a polyhedron with 2A^ — 2 faces, of which 12 are pentagons and the remainder are hexagons. 
Such polyhedra are familiar from the study of fuUerenes in carbon chemistry, with the charge 
17 monopole providing the most famous example of the truncated icosahedron, associated 
with the buckyball. The monopole energy density isosurfaces presented in Figure [2] are 
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qualitatively the same as those for Skyrmions, where the connection with fullerenes was first 
observed [TJ]. 




Figure 2: Energy density isosurfaces for monopoles with charges = 2, 3, 4, 7, 12, 17. Each 
of these monopoles has either a continuous or Platonic symmetry and the plots are displayed 
to scale. 

Note that the difference between the approximate rational map energy (the last expression 
in (12.161) ) and the magnetic bag energy (I2.18P is a factor a/3/2 = 1.22... This suggests that 
generically the energy of the rational map approximation may be as much as 20% above 
the true value. This is much larger than the error in the rational map approximation to 
Skyrmions [5], which is typically a few percent. This means that some caution must be 
exercised when concluding that the symmetry predicted by the rational map approximation is 
indeed the symmetry of the minimal energy monopole. However, the fact that the numerical 
solutions, which share the same symmetries as the rational maps, have energies that tend 
towards the magnetic bag values, suggests that these solutions are strong candidates for 
minimal energy monopoles. 

Further evidence to support this claim is provided by considering alternative rational 
maps. For example, if the axially symmetric degree three rational map R = is used to 
provide the initial (and boundary) conditions for the simulation then an axially symmetric 
3-monopole is obtained with an energy E = 2n x 5.64; this is a little larger than the energy 
E = 27: X 5.54 of the tetrahedrally symmetric 3-monopole. Furthermore, symmetry breaking 
perturbations of axially symmetric maps have been simulated for = 3 and = 4 by using 
the rational map R = + j^z'^~^. In both cases the final numerical solution obtained is a 
Platonic monopole (tetrahedral for = 3 and cubic for A^ = 4), making it extremely likely 
that these are the minimal energy monopoles. It is difficult to perform similar simulations 
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using asymmetric maps for larger values of A^, both because of the dramatically increased 
simulation time required in comparison to symmetric maps, and the fact that the perturbed 
map must remain as the boundary map throughout the simulation, which leads to increased 
angular derivatives as the charge increases. 

Other symmetric maps that are not X minimizing have also been used in simulations. 
For example, there is a degree five map with octahedral symmetry [3] that yields a numerical 
solution for an octahedron with energy E = 2Trx 10.98; this is slightly larger than the energy 
E = 27r X 10.94 of the less symmetric 5-monopole with only symmetry. There is also 
an icosahedrally symmetric degree eleven map [3J that produces an icosahedron with an 
energy E = 27r x 32.85, again larger than the less symmetric D^h 11-monopole with energy 
E = 2tt X 32.49. It seems likely that these additional solutions are saddle points. 

The rational map approximation will predict the correct symmetry of the minimal energy 
monopole providing the error in the approximation is similar for competing local minima. 
One situation in which this is clearly not the case is the flat space limit L — )■ oo of Minkowski 
spacetime, with its 4iV- dimensional BPS moduli space Mn- The rational map approximation 
(12. 6p survives the limit to Minkowski spacetime, but it does not capture the BPS moduli 
space. In this case, the rational map dependence of the energy is not a physical property 
but rather an indication of the failure of the rational map approximation to accurately 
describe a particular monopole in the moduli space. The rational map approximation is at 
its most accurate for shell-like configurations, hence X may be thought of as a measure of the 
deviation from spherical symmetry of an A^- monopole, which is unattainable for > 1. For 
example, for = 4 the X minimizing map has cubic symmetry and the fields (12. 6 p provide 
an approximation to the exact cubic 4- monopole [3], which may be regarded as the point 
in A^4 which is the closest to a spherical 4-monopole. In this example the energy of the 
rational map approximation turns out to be 10% above the BPS energy, and the error is 
obviously greater than this in attempting to describe any other monopole in the moduli space 
Ad^. Despite these limitations, the rational map approximation may turn out to be useful 
in Minkowski spacetime, as it provides explicit monopole fields that can be used in analytic 
approximations, or to provide initial conditions in any numerical computations involving 
monopoles. 

The magnetic bag approximation assumes that the Higgs field vanishes throughout a 
ball. However, exact monopole solutions have a finite number of zeros of the Higgs field. In 
the rational map approximation it is assumed that there is a single zero at the origin, with 
multiplicity A^ for an A^-monopole. This is a correct description of almost all the numerical 
monopole solutions presented in this paper. The exceptions are the monopoles of charges 
three, five and eleven, describing a tetrahedron, octahedron and icosahedron respectively. 
For these monopoles, of which only the A^ = 3 example is minimal energy, the associated 
polyhedron has triangular faces and A^ + 1 vertices. There are A^ + 1 zeros of the Higgs field 
associated with these vertices (but at a smaller distance from the origin) plus an additional 
zero of the Higgs field at the origin with a negative multiplicity (termed an anti-zero). 
This mirrors the situation in Minkowski spacetime, where the same phenomenon has been 
discovered for the tetrahedral A^ = 3 and octahedral A^ = 5 examples [TTJ [18] using the 
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Nahm transform. The = 11 icosahedral example has not been investigated in Minkowski 
spacetime because the Nahm data is not known in this case. 

In Figure [3] the properties of the icosahedral 11-monopole in AdS are displayed. On the 
left is an energy density isosurface, confirming that the energy density is maximal on the 
vertices of an icosahedron. On the right is a plot to highlight the positions of the zeros of the 
Higgs field, by displaying an isosurface where |$| is small; in this particular case |$| = 0.06. 
This plot reveals that there are 12 zeros of the Higgs field on the vertices of an icosahedron 
and an anti-zero at the origin. The surface around the anti-zero appears much larger than 
that around the 12 zeros, which may indicate that the variation of the length of the Higgs 
field is reduced around the anti-zero. However, this conclusion is not certain because of the 
visual distortion in size associated with the AdS metric. 




Figure 3: On the left is an energy density isosurface for the icosahedral 11-monopole. On 
the right (to scale) is the associated isosurface where |$| = 0.06, indicating the positions of 
the zeros of the Higgs field. There are 12 zeros on the vertices of an icosahedron and an 
anti-zero at the origin. 

It has been suggested, at least in Minkowski spacetime, that there may be different types 
of monopole bags depending upon whether the zeros of the Higgs field are located at the 
origin or on the shell of the bag [19] . The numerical results in AdS show that there are indeed 
two different types of monopole solutions characterized either by A^ zeros at the origin or 
A^ -|- 1 zeros near the points of maximal energy density, plus an anti-zero at the origin. In 
AdS it appears that for A^ > 3 the former type of monopole solution has lower energy. 

4 Monopole walls 

In the large A^ limit, monopoles in AdS can be studied by zooming to the Poincare patch 
of AdS. This results in the emergence of a monopole wall, associated with the shell of the 
A^-monopole on which the energy density is maximal. Monopole walls were initially studied 
by Ward [S] in Minkowski spacetime, as novel domain walls in which the magnetic field 
along a line perpendicular to the wall tends to zero on one side of the wall and to a non-zero 
constant on the other side of the wall. In contrast to standard domain walls, the fields are 
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not independent of the coordinates perpendicular to the wall, but are periodic with a spatial 
variation of the energy density and magnetic field. 

In Minkowski spacetime monopole walls have infinite energy per unit area. However, in 
AdS the energy per unit area of a monopole wall is finite, and can be investigated using a 
magnetic bag style approximation [1]. The magnetic bag approximation ignores the spatial 
distribution of the fields within the wall, and is therefore unable to make any predictions 
about the symmetry or lattice structure of the wall. Numerical results in Minkowski space- 
time [H] are unable to shed any light on this issue, as the symmetry is controlled by free 
parameters in the solution, that are essentially inherited from the BPS moduli space for 
finite charge. 

A variant of the rational map approximation is introduced in this section, that is suitable 
for studying monopole walls in AdS. This approximation is refined enough to address the 
spatial distribution of the fields and suggests that the minimal energy per unit area is ob- 
tained from a monopole wall with a hexagonal lattice. Numerical simulations are performed 
that support this conclusion, which is again in agreement with the Skyrme model, where a 
similar hexagonal lattice exists [9]. 

In most applications of holographic methods, the Poincare patch of AdS is the correct 
arena in which to apply the AdS / CFT correspondence. The results of an analysis of monopole 
walls in AdS is therefore likely to be of some interest within the context of holography, and 
in particular for applications to condensed matter systems where magnetic fields play a 
prominent role. 

In the Poincare patch the planar metric of AdS reads 

ds'' = ^(-dt^ + dx^ + dy^) + ^dr^, (4.1) 

where r > is the radial variable in the bulk, with three-dimensional Minkowski spacetime 
obtained at the ultra-violet boundary r — oo. 

Working in a fundamental torus, T^, in the {x,y) plane, the static energy in the volume 

X M is given by 

(4.2) 

The wall version of the rational map approximation has the same form as (12. 6p 

$ = iHn ■ T, Aj = ^{1- K){nx dju) ■ r, (4.3) 

where H{r), K{r) are profile functions and n is again related to a Riemann sphere coordinate 
R through (12.51) . However, in this case R{z) is a periodic function in T^, where z = x + iy 
is the complex coordinate in the plane. 

The boundary conditions on the profile functions are again determined by regularity 
and finite energy to be H{0) = 0, H{oo) = 1, K{0) = 1, K{oo) = 0, where the vacuum 
expectation value |$| — )■ 1 as r — t- oo has been imposed. 
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The ansatz (14 .Sp implies that = and 

F,y = 2t{l-K^)Jn-T, (4.4) 

where 

J=,_ L.,.^'. (4.5) 



1 + |i?P)2 



dz 



The monopole charge in the fundamental torus is the degree of the map from the torus to 
the sphere and is given by 

1 



N = - Jdxdy. (4.6) 
Using the approximate fields (14.31) and performing the integration over yields 

(4.7) 

where A is the area of the torus and 

A f 

X = ^ / J^dxdy, (4.8) 

is a quantity that is independent of the area A. 

As r — )■ 00 the abelian magnetic field perpendicular to the wall is 

B = -lTr(F,,$) = 2J, (4.9) 



with magnetic flux per unit area 



B. = ^l^^2Jdxdy=^-^. (4.10) 

There is a one-parameter family of magnetic walls, labelled by 5^, the magnetic flux per 
unit area at the ultra-violet boundary. Using (14. 7p and the definition (l4.1Up the energy per 
unit area of the wall may be written as 

Before presenting an analysis of the energy (14. lip it is first worth discussing the wall analogue 

of the magnetic bag approximation [1]. 

To the infra-red side of the wall r < a, the profile functions are taken to be = and 
= 1, so that the energy density vanishes. To the ultra-violet side of the wall r > a, the 

simplification is to assume that K = 0, so that the energy per unit area becomes 



E rrrrAh^„^i|i/Ly|/.y^^_ ^^^^^^ 



r 

A {\Lj 2" ' 2m\r J j \L 
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where once again the contribution on the wall r = a has been ignored. 

The energy minimizing profile function, satisfying the boundary condition H{(t) = 0, is 
easily found to be 



H = l-—. (4.13) 



The position of the wall is 



with a resulting energy per unit area 



3N 



E /4S3LX3/2 



A V 3A^3 



(4.14) 



(4.15) 



The final approximation required to reproduce the result of Bolognesi and Tong [Ij is to 
assume that the magnetic field is independent of the coordinates parallel to the wall. This 
simplification implies that B = B^, = 2J is a constant, and hence the approximation X = N"^. 
This yields the wall position and energy per unit area as 



Returning to the energy (14. lip , an explicit expression for R{^z\ the map from the torus to 
the sphere, is required to make further progress. To obtain a wall with hexagonal symmetry 
this map is chosen to be proportional to the Weierstrass elliptic function ip{z) defined by the 
equation 

p'2 = V_4. (4.17) 

This elliptic function has periods Wi = r(|)r(|)/(2-\/37r) and iiJ2 = iiJ\e™l'^^ giving the 
required 60° angle between the fundamental periods. The precise form taken for the map is 
Ri^z) = cp{z/a), where a and c are real constants, with a determined in terms of the area A 
of the torus by = 2A/{\/3ujf). Recall that I/N"^ is independent of A (and hence a). 

The elliptic function has a double pole in its fundamental parallelogram and describes 
a map with degree N = 2. Integrating over the torus reveals that I/N"^ is minimized for 
c = 0.70, when it takes the value X/N"^ = 1.21. Note that this value is consistent with the 
large limit of the minimizing rational maps discussed earlier. Using the value X/iV^ = 1.21 
allows the energy per unit area (14. lip of the elliptic map approximation to be calculated by 
computing the minimizing profile function. 

There is a scaling isometry of the metric (14. ip that relates monopole wall solutions with 
different values of the magnetic flux per unit area. Explicitly, the metric is invariant under 
the scaling x i— )■ X~^x, y t— )■ X~^y, t i— >■ X~^t, r ^ \r which yields B^ i— )■ X^B^ and 
E/A I— )■ \^E/A. This scaling symmetry may be used to restrict the computations to a 
convenient positive value of S*. 

Numerical field theory simulations allow a computation of the monopole wall energy (14.20 
using a simulated annealing algorithm similar to that discussed in the previous section. It 
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is convenient to perform the simulations over two copies of the torus, so that a rectangular 
domain {x,y) G [0,au;i] x [0, \/3aui] may be used with periodic boundary conditions. The 
numerical grid contained 52 x 90 grid points to cover each rectangle. 

It is numerically more efficient if all the spatial coordinates take values in a finite range, 
hence the simulations use the variable u = r/(l + r), taking values in the unit interval. This 
interval is covered with 50 grid points, so the total grid contains 52 x 90 x 50 points. The 
initial conditions are taken from the elliptic map approximation. 

In Figure m the image on the left is an energy density isosurface for the numerical solution 
with = 10, where the hexagonal structure of the lattice is clearly visible. As the region 
contains two copies of the torus then the charge in the displayed region is iV = 4. 




Figure 4: Energy density isosurfaces for monopole walls with B^, = 10. On the left is the 
hexagonal wall and on the right is the square wall, which has a slightly higher energy. 

The energy per unit area of this hexagonal wall is computed to be E/A = 39.1. For 
comparison, the energy per unit area of the elliptic map approximation is E/A = 42.7 and 
for the magnetic bag style approximation is E/A = 36.5. As in the case of monopoles in 
global AdS, this demonstrates that the magnetic bag style approximation provides a lower 
bound on the energy, to complement the upper bound of the elliptic map approximation. 

In Figure [5] the length of the Higgs field |$| is plotted as a function of r/(l + r) for 
the case -B^, = 1. The dashed curve is the elliptic map approximation and the dotted curve 
is the magnetic bag style approximation. The solid curve is the result from field theory 
simulations, plotted along a generic line perpendicular to the wall. It can be seen that the 
elliptic map approximation is in excellent agreement with the field theory simulations, and 
the magnetic bag approximation provides a good description everywhere except around the 
actual wall itself. 

A monopole wall with a square lattice can be obtained by using the Weierstrass elliptic 
function p{z) defined by the equation 

p'^ = Ap^-Ap, (4.18) 

which has periods Ui = T{j)'^ /{2^/2tt) and UJ2 = ioJi-, producing a 90° angle between the 
fundamental periods. 
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Figure 5: |$| as a function of r/(l + r). The dashed curve is the elhptic map approximation, 
the dotted curve is the magnetic bag style approximation, and the sohd curve is the result 
from field theory simulations. 



Taking a map of the form R{z) = cp{z/a), with = A/uf, reveals that I/N"^ is min- 
imized for c = 1.00, when it takes the value X/N"^ = 1.30. This is greater than that of the 
hexagonal lattice, and supports the view that a hexagonal architecture produces minimal 
energy. 

Further evidence is provided by using the approximate fields of the square wall as initial 
conditions in the field theory simulations. To facilitate a comparison with the simulations of 
the hexagonal wall, two copies of the torus are again taken, by using the rectangular domain 
{x,y) G [0,aa;i] x [0,2aa;i], covered by 50 x 100 grid points. In Figure H] the image on the 
right is an energy density isosurface for the resulting numerical solution with B^, = 10. The 
energy per unit area of this square wall is E/A = 39.5, which is slightly greater than the 
value E/A = 39.1 for the corresponding hexagonal wall. 

The connections between monopoles and Skyrmions suggests a qualitative understanding 
of the preference for a hexagonal monopole wall, since it is known that this is the minimal 
energy form for a wall in the Skyrme model [9]. It seems that the same generic structure 
arises when both types of soliton overlap, even though the source of the soliton attraction 
is different in the two theories; there being an attractive force between two solitons in the 
Skyrme model and an attractive force between a soliton and the origin of AdS in the case 
of monopoles. 




18 



5 Conclusion 



Motivated by applications in holographic theories, a detailed analytic and numerical study 
has been performed for monopoles in AdS, and the associated monopole walls that arise in 
the large charge limit. An approximation has been introduced, using rational maps between 
Riemann spheres, and it has been confirmed that this provides a reasonable description of 
the fields when compared with the results from field theory simulations. The results of these 
simulations also confirm that the magnetic bag approximation yields energies that are close 
to the true monopole energies, even for reasonably small values of the charge. In all the 
numerical computations presented in this paper the value L = 1 has been used as a generic 
choice for the AdS radius. It might be interesting to extend the computations to other values 
of L, or equivalently to a range of vacuum expectation values of the Higgs field. 

The arguments presented in [T] suggest that, for certain regions of parameter space, 
the monopole wall (with its lattice structure) may be favoured over the more conventional 
magnetically charged black hole, which is homogeneous at the ultra-violet boundary. It 
would be interesting to investigate this issue by extending the methods described in this 
paper, for a fixed AdS background, to the situation including gravitational backreaction and 
finite temperature. 

Finally, the connection between monopoles and Skyrmions in Minkowski spacetime has 
been enhanced by considering monopoles in AdS. A natural extension would therefore be to 
study Skyrmions in AdS, to make a comparison with the results presented here for monopoles. 
The rational map approximation easily extends to Skyrmions in AdS and yields the same 
functional on the space of rational maps. This would seem to imply that monopoles and 
Skyrmions in AdS have a similar form. However, there is a caveat to this conclusion. It 
appears that, as far as Skyrmions are concerned, the curvature of hyperbolic space plays a 
similar role to that of a pion mass [20]. It is known that Skyrmions with massive pions, in 
Minkowski spacetime, are shell-like (and hence described by the rational map approximation) 
only for sufficiently low baryon numbers, and take a different form, including clusters, above 
a critical value [21] 122] . It may therefore be possible that Skyrmions in AdS could take more 
exotic forms than monopoles, such as clusters or multiple shells. It might be interesting to 
explore these possibilities. 
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